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In this work we re-examine the correlation between B — mixing and Lepton Flavour 
Violation in the light of recent experimental measurements in the B s system. We perform a generic 
SUSY analysis of the allowed down squark mass insertion parameter space. In the SUSY GUT 
scenario this parameter space is then used to make predictions for LFV branching ratios. We find 
that the recent measurement for the CP phase <j) s excludes the lowest rates for r — > /17 and provides 
a lower bound of ~ 3 x 10 -9 for tan f3 — 10. Future experimental improvements in the bound on 
t — *• y/y and the measurement of <j> a will constitute a strong test of the SUSY GUT scenario. 



PACS: 12.10.-g; 12.60. Jv; 13.35.Dx; 14.40.Nd; 13. 20. He 



Nucl.Phys.B802:63-76,2008 



I. INTRODUCTION 

In recent years there have been great efforts made in 
the search for New Physics (NP) in the B system which 
has been seen as an extremely promising place to find NP 
effects as deviations of flavour-changing neutral-currents 
(FCNC) from their Standard Model(SM) expectations. 
The first measurements of B s — B s mixing have been ob- 
served at D0 [TJ and CDF [3] experiments reporting a 
mass difference, 

17 ps" 1 < AM S < 21 ps" 1 (90%C.L. D0) (1) 

AM, = 17.77 ± 0.10 ± 0.07 ps -1 (CDF) (2) 

The constraints from eq. Q strongly restrict the allowed 
parameter space for the NP contribution to B s —B s mix- 
ing ei ei hdi nn us na na nsi - 

The first measurement of the CP phase associated with 
B s — B s mixing, tj) a , was reported by the D0 collabo- 
ration [Sj. More recently, flavour-tagged measurements 
have now been made by both the CDF collaboration [4 
and the D0 collaboration [5]. Combining these latest 
results the UTfit group 6 found that the new physics 
phase deviates from zero by about 3cr, 



S = (<f>» p - 20s) 



(3) 



with p s = 0.018 ± 0.001 and, 

p = [-60.9, -18.58]° U [-156.90, -106.40]° (4) 

at 95% C.L. There still remains a two fold ambiguity in 
the measurement due to the symmetry in cj) s — * (it — cp s ). 

In the Bd system the measured mass difference given 
by the "Heavy Flavour Averaging Group" (HFAG) [17 is, 



AM d = 0.507 ± 0.004 ps" 1 



(5) 
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and the phase, 0<j, has been measured as, 

(j) d = 0.757 ± 0.044 (6) 

A very appealing model of NP is the supersymmetric 
grand unified theory (SUSY GUT). In such models there 
is a deep underlying connection between quarks and lep- 
tons above the GUT scale. For example in SU(5) the 
singlet down quarks are related to the lepton doublet as 
they live in the same 5 representation. This deep relation 
between down quarks and leptons is broken at the GUT 
scale but its presence may still be felt all the way down 
at the Electroweak scale. As a result of this relationship 
the right-handed down squark and left-handed slepton 
mixings are related to each other. In such SUSY GUTs 
there is hence a deep connection between LFV rates such 
as BR(r — * /ry) and FCNCs such as B s —B s mixings and 
its associated CP phase [18 l [T9 l l20 l 121] . 

In this work we re-examine the correlation between FC- 
NCs and LFV rates in SUSY GUTs. Specifically we are 
interested in the connection between the allowed NP pa- 
rameter space of B s —B s mixing and the rate of t — * (j/y- 
To this end we first perform a generic SUSY analysis, of 
the allowed parameter space of down squark mass inser- 
tions <5^jj. Our analysis in generic in that no particular 
SUSY-breaking mechanism shall be assumed. Having de- 
termined the allowed parameter space we then exploit the 
GUT relationship of squark and slepton mass insertions 
to make realistic predictions for the rates of r — > (X'j, 
t — * e7 and their ratio. Recently a similar study was 
undertaken [TO] on the correlation between LFV and B q 
mixing. We believe that the NP contribution to Bd mix- 
ing was under-estimated in that work leading to over- 
estimates for the branching ratio of r — > ej. In this work 
we use an analytic form for the gluino contribution to 
Bd mixing derived in appendix A in the mass insertion 
approximation. This allows us to correctly compute the 
NP contribution to both Bd and B s mixing. Further- 
more the recent measurement of the CP phase <j) s has 
prompted a re-examination of this correlation. Crucially 
we find that these recent measurements severely restrict 
the allowed parameter space and provide a lower bound 



2 



for the rate of r — > fij. Finally, we explore an exam- 
ple of an SO(10) inspired model and discuss the tension 
between the constraints of LFV and FCNCs. 



II. CONSTRAINTS ON THE SUSY 
CONTRIBUTION TO B„ MIXING 



The AB = 2 transition between B q and B q mesons is 



defined as, 



(B° q \nt f B f = 2 \B° q ) 



2M Bg M? 2 



(7) 



where M Bq is the mass of the B q meson. We can then 
define the B q mass cigenstate difference as, 



AM q = M q H -M q L = 2\M? 2 
and its associated CP phase, 

0, = arg{Ml 2 ) 
In the Standard Model M\ 2 is given by, 



(9) 



Mff M = ^fM Bq if f Bq B Bq {V; q V tb f S (x t ) (10) 

where Gp is Fermi's constant, M\y the mass of the 
W boson, fj B = 0.551 is a short-distance QCD correc- 
tion identical for both the B s and Bd systems. The 
bag parameter B Bq and decay constant f Bq are non- 
perturbative quantities and contain the majority of the 
theoretical uncertainty. V tq and V t b are elements of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and 
S (x t = fh 2 /M^) = 2.32 ± 0.04. 

The NP contribution to B q mixing may be parameter- 
ized in a model independent way as, 



AM, = 



AA/ SM h + R 



6 sm 



6 SM 



arg(l 



(11) 
*<) (12) 



where R a 



r a e 



■ Mf 2 NP /Ml 2 M denotes the relative 
size of the NP contribution. The dominant SUSY contri- 
bution to B q mixing comes from the gluino contribution 
which may be written as, 



R\ = al(m- g ,x) [(6 d q3 ) 2 RR + (6 d 3 ) 2 LL ] 



+al(m g ,x)(S d q3 ) LL (6 d q3 ) RR + ... (13) 

Here we have ignored terms proportional to S RL LR mass 
insertions as they are expected to be heavily suppressed 
due to constraints from b — > [18 . In appendix A we 
give details of the functions a\ and a\. 

We can now constrain both the magnitude and phase of 
the NP contribution, r q and a q , through the comparison 
of the experimental measurements with SM expectations. 



From the definition of eq. (Ill we have the constraint 



Pi 



AM q 

AM| M 



1 + 2r q cos a q + r 2 



(14) 



for r q and a q . The values for p q given by the UTfit 
analysis [6, 7 at the 95% C.L. are, 



Pd 

Pa 



[0.53, 2.05] 
[0.62, 1.93] 



(15) 
(16) 



From eq. (12) the phase associated with NP can also 



be written in terms of r q and <j q , 



/NP 



COS 



/NP 



1 + 2r q cos <T q + r 2 q 

1 + T q COS (Tq 

1 + 2r q cos o q + r 2 



(17) 



(8) Here [SIE] gives the 95% C.L. constraints, 



^ p = [-16.6, 3.2]° 

iNP 



= [-60.9, -18.58]° U [-156.90, 



(18) 

-106.40]° (19) 



In order to consistently apply the above constraints 
all input parameters arc chosen to match those used 
in the analysis of the UTfit group [51 [7] with the non- 
perturbative parameters, 



f Bs B Bs = 262±35MeV 

I = 1.23 ± 0.06 MeV 

= 230 ±30 MeV 

= 189 ±27 MeV 



1b. 
1s d 



(20) 
(21) 
(22) 
(23) 



Now we have a clear picture of the constraints to be 
imposed on the r q — o q plane as listed in eq. ( 14p9 1. In 
the following sections these constraints shall be used in 
a generic SUSY analysis of the mass insertion parameter 
space, followed by a study of the correlation of B q mixing 
and LFV rates in the minimal SUSY GUT. 



III. CORRELATION BETWEEN B q MIXING 
AND LFV RATES IN SUSY GUTS 

In many ways the prototypical GUT is SU(5) as it is 
the minimal group which can unify the gauge group of 
the Standard Model. In the SU(5) GUT the quarks and 
leptons are placed into 10 = (Q, u c , e c ) and 5 = (L, d c ) 
representations. Due to the symmetry of this simple 
SUSY GUT there exists relations amongst the slepton 
and squark soft SUSY breaking masses, 



TO I> 



(24) 



These relations hold for scales at and above the GUT 
scale. Interestingly this implies that left-handed slepton 
mixing and right-handed down squark mixing are related. 
This relation can still be felt very strongly at the Elec- 
troweak scale in the correlation of LFV rates and FCNCs. 

Further renormalization group(RG) evolution down to 
the Electroweak scale has very little effect on the size 
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of these off-diagonal elements [13]. So we may assume 
that these GUT scale values are approximately equal to 
their values at the clcctrowcak scale. Hence it is a fair 
approximation to assume that (to 2 ^)^- ~ (m^)ij at the 
electroweak scale. Then the (8fj)n.R = kp^)%jl m \ con- 
tributions to FCNC and (fiy)££ = (m|)y/m| contribu- 
tions to LFV are clearly correlated. We may explicitly 
write the rate of Zj — > L7 in the very suggestive form [19j . 



Bi?(Z, -> l 3l ) 



cr m 



G 2 F Mg 



/3 



(25) 



where m§ is the average squark mass and Ms is the typ- 
ical SUSY scale. 

We also need to consider the RGE effects of the CKM 
mixings in the left-handed down squark matrices. These 
effects can be approximated as, 



(4) 



~ 1 v 2 v*v ( 3m o 



Al) 



to? 



In 



M* 



M- 



(26) 



Here Too is the universal scalar mass, Aq the universal A- 
term and M* is the scale at which the flavour blind soft 
SUSY breaking is communicated. These effects may be 
quite important in the SllSrr contribution to eq. (13 1 
due to the typically large value of a 4 ss — 100 Oi- We 
shall take a minimal approach and assume M* to be the 
GUT scale ~ 2 x 10 14 GeV, in which case the mass in- 
sertions are of the order, (5 d 3 )LL ~ A 2 and (S d 3 )LL ~ A 3 . 
Larger values of M* will lead to an even more restricted 
parameter space for (Sf 3 23) RR- 



If we assume that the second term of eq. (13 1 domi- 
nates then we may derive the relation, 



(<&) 



LL 



X 



(6 d 13 )LL 

2 BR[T-^nj) 
BR(t -> ei) 



RR 



(^13) RR 



(27) 



here we have made use of the ratio, 1 0.4 1 / la.4 1 ~ Vtd/^ts ~ 
A 2 , as derived in appendix A. This relation differs from 
that derived in [19] where it was assumed that, a\ja% = 1. 
In appendix A we show the full form of the functions 
di and 04 where it is shown that the correct ratio is 
l a ll/l a 4l ~ Vtd/Vts ~ A 2 - As a result the allowed param- 
eter space for the mass insertion (Si 3 )rr is suppressed 
relative to (5$ 3 )rr by a factor A 2 as shown above. 



IV. NUMERICAL RESULTS AND DISCUSSION 

In this section we shall show the results of numerical 
calculations of our generic SUSY analysis of the allowed 
parameter space for the mass insertions (Sf 3 2 3)RR- We 
then make use of this allowed parameter space to study 
the correlation of FCNCs and LFV rates. 

For the numerical analysis we fixed the values of 
tan/3 = 10, M s = m Q = M 1/2 = mg = (300, 500) GeV, 



Aq = and the ratio x = nig/m- — 1. Taking these val- 
ues we scan over (^3 13) rr and require fits to the values 
of p q and (p^ F as described in the previous section. 
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FIG. 1: Allowed parameter space for the mass insertion 
(S23) rr, with m s = 300(blue/pink) and 500(red/green) GeV 
respectively. Red/Blue points are constrained by the mea- 
surement of AM S while green/pink points have the extra con- 
straint from the measurement of <j!> s both at the 95% C.L. Solid 
black lines show the respective constraints for each value of 
mo from bounds on r — > fi'y from Belle [24] . 

In fig. [T] we show the allowed parameter space for the 
mass insertion (<5;>3)-R-R as dictated by the constraints im- 
posed from AM, (red/blue) and <j) s (green/pink). The al- 
lowed regions form rings in the complex ($23) RR plane 
with wedges cut into it corresponding to allowed regions 
from the measurement of <fi s . It is clear therefore that 
the recent measurements of (f> s substantially restrict the 
allowed parameter space and provide important informa- 
tion. In the SUSY-GUT scenario the same mass insertion 
also contributes to the LFV decay r — > /ry, the resulting 
bounds from which can be seen by the solid black line 
shown in fig [T] The two black lines correspond to the 
values too — 300, 500 GeV and clearly exclude consider- 
able regions of the (^23)^^ parameter space. In fact one 
of the two allowed regions is strongly disfavoured in the 
SUSY-GUT scenario. 

Fig. [2] shows the respective allowed parameter space of 
the mass insertion (Sf 3 )RR constrained by the measure- 
ments of AMrf and <t>d- For (Sf 3 )RR we can see that the 
measurement of 4>d has reduced the allowed parameter 
space to a small slice in the complex (Sf 3 )RR plane. In- 
creasing the size of rriq leads to an increase in the allowed 
parameter space due to a suppression of the functions 
ai, 4 ~ l/m|. 

From these allowed regions of parameter space we can 
plot the corresponding allowed regions of the r q — a q 
parameter space. These are shown in fig. [3] and fig. [4] 
Without the restrictions from the CP phases 4>d,s the al- 
lowed regions of parameter space form loops in the r q — a q 
plane. The constraints from the CP phases represent 
slices as shown in [5) Q2] . The accurate measurement of 
4>d leaves the allowed region of parameter space as a thin 
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slice shown by the green points in fig. [3] On the other 
hand, the present measurement of 4> s has a two fold am- 
biguity shown by the two green wedges in fig. [4j 
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FIG. 2: Allowed parameter space for the mass insertion 
13) RRi with rrig — 300(small) and 500(large) GeV respec- 
tively. Red/Blue points are constrained by the measurement 
of AMd while green/pink points have the extra constraint 
from the measurement of <fid both at the 95% C.L.. 

Next we would like to discuss the implications of this 
allowed parameter space on predictions for the LFV rates 
of r — > (j/y and r — * ery. We have shown that in a SUSY 
GUT the mass insertions associated with the off-diagonal 
down squark mass matrix elements and those of the slep- 
tons are clearly related. This relation leads to a strong 
correlation between the B q mixing and LFV rates. Here 
we would like to expose this correlation in making realis- 
tic predictions for LFV rates and to extract information 
about the NP CP phase. 
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FIG. 3: Allowed region of Vd — ad parameter space for the NP 
contribution to Bd mixing. Red points agree with the AMd 
constraint while green points also agree with that of cj>d both 
at the 95% C.L.. 

Let us first consider the 23 sector. Fig. [5] shows the 
predicted rate for r — > fij plotted against the NP phase 
<fi^ p . The possibility of 4^ p = is excluded at the 3er 



level [BJ, which also excludes the lowest rates for r — * fvy 
at the centre of the plot where </>^ p — > 0. The gen- 
eral feature of the plot is that larger LFV rates are pre- 
dicted for larger NP phases. Without the input of <j> s 
this plot is symmetric in (p^ p — ► — (f)f and extends for 
all <fi^ p € (— 7r, 7r). The inclusion of the 4> s constraint 
drastically reduces the allowed parameter space and the 
resulting prediction space, as shown in fig. [5| These pre- 
dictions are independent of the choice of Ms = rrig due 
to the dependence of both BR(Zj — > Ijj) and the 

functions a\ t i shown in appendix A. There are two al- 
lowed regions, one with large NP phase and one with 
smaller NP phase. From fig. [5] it is clear that the up- 
per bound on r — > /17 disfavours the large phase solu- 
tion. For the second region, there is a lower bound on 
the BR(r — ► hj) > 3 x 10 -9 . This second allowed re- 
gion is just below the present experimental bounds for 
the branching ratio, BR(t — > /17) < 6.8 x 10 -8 BaBar 
(2SJ|S3] and BR(r -> ^7) < 4.5 x 1CT 8 Belle [M]. From 
these plots it is clear that improvements of the bound on 
t — ► pr/ in conjunction with an improved measurement 
of (f> s will provide a much stricter test of the SUSY GUT 
scenario. 
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FIG. 4: The allowed region of the r s — a s parameter space 
for the NP contribution to B s mixing. Points in Red agree 
with the constraint from AM S , while green points also agree 
with the 4>s measurement both at the 95% C.L.. 

In the 13 sector, we can see from fig [2] that the al- 
lowed parameter space for the mass insertion ((5f 3 )_Ri{ is 
far more restricted. This is due to the smaller mass differ- 
ence in the Bd system, as a result the functions a\ and 
are enhanced by the ratio, AM d /AM s » V t 2 d /V t 2 s ~ A 2 , 
see appendix A. This restriction leads to a much more 
suppressed prediction for the rates of r — > 67. These 
rates are plotted in fig. [6] against the allowed values of 
the NP phase </>^ p . Including the constraints from <f>d 
the predicted values are in the region of 10~ 14 — 1CP 9 . 
The present experimental bounds for the branching ra- 
tio are BR(r -> e 7 ) < 1.1 x 1CT 7 BaBar [22| |23] and 
BR(r -> e 7 ) < 1.2 x 1(T 7 Belle (24]. Looking at fig. § we 
can see that the allowed region of the (5f 3 )_R_R parameter 
space predicts rates well below the present experimental 
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bounds. Therefore the present bounds have no impact 
on the allowed parameter space. 




FIG. 5: Predictions for BR(r — > fij) from constraints on B s 
mixing. The red points conform to the AM S constraint while 
green points show the additional constraint from <j> s , both at 
95% C.L. The present experimental bound is shown by the 
solid black horizontal line. 



The ratio of the branching fraction of r — * /17 and 
t — ► c^/ is less dependent on the SUSY parameter space 
and so it is interesting to consider the allowed size of 
this ratio. Using the parameter space of (8f 3 )jni and 
(S$ 3 )rr allowed by the constraints of pd, s and (f>^ we can 
make predictions for the ratio BR(r — > [i-f) /BR(r — > ej). 
Fig. [7] shows the resulting plot of BR(r — > /i7)/BR(Y — > 
e-f) plotted against the prediction for the NP CP phase 
c^ p . We can see that the numerical results show this ra- 
tio ranging from 0.01 up to tens of thousands. Applying 
the constraints from the CP phases <pd,s implies that this 
ratio must lie in the region > 10. This large ratio is again 
due to the suppression of the (Sf 3 )uR allowed parameter 
space due to the smallness of the mass difference in the 
Bd system as opposed to the larger mass difference in 
the B s system. Here our results disagree with those pre- 
sented in [19] due to their erroneous assumption of the 
ratio af/af = 1. 

Let us now look at a specific example of a simple 
SUSY GUT. The example we pick is an SO(10) inspired 
model as introduced in [201 ■ In this model there are 
two Higgs representations one for the up/neutrino and 
down/charged lepton Yukawa couplings separately. The 
neutrino Yukawa coupling is then related to the up quark 
Yakawa, where two limiting cases are possible; (i) Y v has 
CKM like mixing, (ii) Y v has MNS like mixing. It is 
assumed that SO(10) breaks to SU(5) at the scale M10 
and further breaks to the SM gauge group at the scale 
Mqut- Assuming that the flavour blind SUSY breaking 
occurs at the scale M* = M10 = 10 17 GeV, then the mass 
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FIG. 6: Predictions for BR(r — > &y) from constraints on Bd 
mixing. 



insertions can be written as, 



(4) 
(*«) 



RR 



LL 



3 2 
■ 8^2 »« Wain 



Al- 



io 



M GU T 

Mio 



Mr 



(28) 
(29) 



where the matrix V is Vckm and Umns f° r case (i) and 
(ii) respectively. The scale Mr ~ 10 15 GeV is the scale of 
the right-handed neutrinos. Case (i) produces relatively 
small LFV rates with (5 l 23 ) LL - 0.008 and (5$ 3 ) RR ~ 
0.003. The present bound on the rate of r — > ^7 leads 
to the weak bound mo > 200 GeV. Alternatively the 
contribution to B s mixing is also very small, for mo > 
200 GeV we have too small a value |i? s | < 0.2. Case 
(ii) produces large LFV rates with (<523)ll ~ 0.25 and 
(^23) RR ~ 0.08. Here the present bound on the rate of 
t — > /i^f leads to the strong bound mo > 1000 GeV. Here 
again the contribution to B s mixing is very small, for 
™o > 1000 GeV we once again have too small a value 
\R,\ < 0.2. 




90 135 180 



FIG. 7: BR(r — + ^7)/BR(r — > cy) from constraints on B s 
and Bd mixing. Red points are from the AMd, a constraint 
only, green points are for both AMd, s and <j>d )S constraints. 
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V. SUMMARY 

In this work we have re-examined the constraints im- 
posed on the parameter space of NP contributions to 
B s and Bd mixing in the light of recent measurements 
of AM S and its associated CP phase <f> s . These con- 
straints were then imposed to make a generic SUSY anal- 
ysis of the allowed parameter space of the mass insertions 
(Sf 3 )RR and (S d 3 ) R R. 

In SUSY GUTs there is a deep underlying connection 
between (s)quarks and (s)leptons. In such a theory we 
should therefore expect that the flavour mixings observed 
in the quark and lepton sectors are correlated. Assuming 
that such a SUSY GUT exists, we have re-examined the 
correlation between LFVs and FCNCs. 

From our numerical analysis we have found that the 
allowed parameter space for the mass insertion (Sf 3 ) RR 
is particularly restricted by the present measurements of 
AMd and cf>d- This is due to the small mass difference in 
the Bd system compared to the B s system. As a result, 
the predicted branching ratio for r — > ej is particularly 
suppressed. On the other hand, the present experimental 
determination of AM S and particularly <fi s are not yet so 
restrictive on the allowed values of (S^rr- The larger- 
mass difference in the B 8 system also means that the 
predicted values of the branching ratio of r — > /17 are 
much larger and close to the present bounds from the 
B factories. The main conclusion of this work is that 
the recent measurement of <j> s substantially restricts the 
allowed parameter space and excludes the region where 
the lowest rates for r — > ^7 are to be found, leading to 
a lower bound of ~ 3 x 1CP 9 for tan/3 = 10. Future 
experimental improvement of the bound on the decay 
rate r — * fij and the measurement of <f> s will lead to a 
strong test of the SUSY-GUT scenario. We have also 
considered a specific example of a SUSY SO(10) GUT 
and found there to be significant tension between the 
constraints from FCNCs and LFV decay rates. 



APPENDIX A: B° A s - B° d s MIXING IN SUSY 
MODELS 



To examine the contributions from new physics we may 
write the AB — 2 effective Hamiltonian in the following 
form, 



We shall also need to define the hadronic matrix elements 
of the above operators such that, 



(B q \Qi\B q ) 
(B q \Q 2 \B q ) 
(B q \Q 3 \B q ) 
(BMB q ) 
(B q \Q 5 \B q ) 



-^MbJIb^) (A7) 

~RB q M Bq fl q B 2 (») (A8) 

--R^MbJIB^) (A9) 

-^RbMbJIB^) (A10) 

-^RbMbJIB^) (All) 



where we have defined Rb = ( — r-vr^—r^ I ■ The val- 
ues of the bag parameters have been calculated on the 
lattice and can be found in [25]. 

The dominant SUSY contribution is from the gluino, 
which gives the following Wilson coefficients [271 1251 

EBEHE2I 



C'l (Ms) = 



216 m? 



24xf 6 (x)+66f 6 (x) (^ 3 ) 
a 2 s (M s ) 



d \2 
LL 



(A12) 



Cf (M S ) = - 



216 to? 



50ixf 6 (x) - 72f e (x)\ (6%) LL (S* 3 ) RR (A13) 
a 2 s (M s ) 



Cl(Ms) = 



216 to| 



24x/ 6 (x) + 120/e (x) (S d q3 )LL(S d q3 ) R R (A14) 



Here only the LL and RR terms have been kept as the 
contributions from RL and LR are suppressed. The 
above Wilson coefficients are evaluated at the scale M5 = 
(M- g + M q )/2 = (V^+ l)M q /2. 

These Wilson coefficients need to be renormalization 
group (RG) evolved to the scale of the bottom quark. 
Using the magic numbers from [32] , we have, 

C M -EE ^ + V cf )/3) ) V ai Cp{M s ) (A15) 



(Al) 



with r\ = a s (Mg) /a s (fit)- Let us now write the B s — B s 
mixing parameter 



The operators are defined as, 



Qi = ^mh^L 
Q2 = <m<?A 

Q 4 = q a R b a L f L b R 
= <%b{(f L b R 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 



AM„ 



2\MU 



2 M 



irq,SU 



= AM? 11 



12 
SM 



(l + R q 



R 9 q 



(A16) 



where R q = r q e 



= Mf^ F /Ml^ M parameterizes the 
contribution from new physics. We consider the domi- 
nant effect to come from the gluino i?| = M\£ /Mfy 
From the above effective Hamiltonian, Wilson coeffi- 
cients, operator matrix elements and RG evolution we 
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can write, 



Rl = al(m S) x) [{5%) 2 RR + {5%)\ 

+al(m 9) x) (S d q3 ) LL (S d q3 ) RR +■■■ (A17) 

have we have only kept the LL, RR terms as they will 
dominate. The complete form of the functions ai(mg, x) 
and ai{m,g,x) are as follows, 



Mf 2 SM a\{mg,x) = 



-ag(M s ) 
216 to? 



A 



\m B JlBMY.^' 1) +^^ 1) )^ i ( A18 ) 



Mf 2 SM a\(mg,x) 



-a 2 s (M Sj ^ „ 2 



216 to? 



m B s f Bs R B q 



i 

i 
i 



(A19) 



Where we have defined, A = 24 x fe(x) + 66 fe(x) and 
B = 504 a: / 6 (a;) - 72f 6 (x) where, 

, , , 6(1 + 3x) lnx + a; 3 - 9a; 2 - 9x + 17 , . nn . 
M*) = J 6 (x-l) 5 (A20) 

~ e(x) Jx(l + x)lnx^-9xi + 9x + l ^ 

and a; = to? /to?. 

The dominant contributions to the above functions 
come from b[ hl) = 0.865, b{ 4A) = 2.87, 6 4 4 ' 5) = 



0.961, &f' 4) 



0.09, and b 



(5,5) 



0.863, with cii 



(0.286, -0.692, 0.787, -1.143, 0.143). 

So we can write the approximate functions as, 



Mff M a1(m g ,x) f 

xA 0.865 ?7 - 286 



a 2 s (M s ) 1 



216 to? 3 



m Bj B Bi(n) 



(A22) 



M 04(7715, x) 



-a 2 s (M S 
216 to? 



[B 2.87 rT 1143 +yl 0.961 7? 

"ijo.ogr?- 1 143 



RB q m Bs fl 3 
- 1143 ] ^ 4 ( M ) 



^ 0.863 77 0143 ] -B B (/x) 



(A23) 



We can see that the values of the functions oy arc not 
the same for the B s and Bd systems. From eq. ( A18|A19 1 
we can approximate the ratio of these functions as, 



\a d A ~ AMf M ~ I/ 2 



(A24) 



so that the value of the functions a\ and 04 are enhanced 
in the Bd relative to those derived in the B s system. This 
enhancement will lead to a more constrained parameter 
space of the 13 mass insertion relative to the 23 sector. 
It is also important to notice that a\ and 04 are complex 
parameters. This is particularly important in the case of 
the Bd system where the phase of M d ^ M , = 2(3, is 
large. 
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